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A second order, O(a^), calculation in perturbative quantum chromo dynamics 
is presented of the longitudinal spin dependence of the cross section for the two 



^ \ particle inclusive reaction p + p ^ ^ + c + X for large values of the transverse 

^ momentum of the prompt photon and charm quark. Differential distribu- 

"X' 

5h ' tions are provided for the spin-averaged cross section and for the two-spin 

longitudinal polarization asymmetry All at the energy of the Brookhaven 
Relativistic Heavy Ion Collider. An assessment is given of the prospects for 
determination of the spin dependence of the charm quark density. 
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I. INTRODUCTION 



Because they couple in a point-like fashion to quarks, the observation of photons with 
large values of transverse momentum in a high energy hadron collision has long been re- 
garded as an incisive probe of short distance dynamics. The dominance of the Compton 
subprocess g + t^— i^7 + gat both leading- and at next-to-leading orders |T]] in perturbative 
quantum chromodynamics [QCD] makes spin-averaged inclusive prompt photon production 
the reaction of choice for investigations of the magnitude and Bjorken x dependence of the 
gluon density of the incident hadrons, G{x, n). The Compton subprocess also dominates the 
dynamics in longitudinally polarized proton-proton reactions as long as the polarized 
gluon density AG{x,fi) is not too small. As a result, two-spin measurements of inclusive 
prompt photon production in polarized pp scattering should constrain the size, sign, and 
Bjorken x dependence of AG{x,fi). 

In previous papers we investigated the associated production of a prompt photon along 
with a charm quark [^Q, p + p^'y + c + X, in spin-averaged reactions at high energy. 
We showed that collider data from this two-particle inclusive reaction at large values of 
transverse momentum may be used to measure the spin-averaged charm quark density in 
the nucleon. In this paper, we extend our previous work by including a full next-to-leading 
order treatment of the longitudinal spin dependence. We specialize top + p^'j + c + Xa.t 
the energies of the Brookhaven Relativistic Heavy Ion Collider (RHIC) facility. 

For values of the transverse momentum p^ of the charm quark much larger than the 
mass rric of the quark, only one direct hard scattering subprocess contributes in leading 
order: the quark-gluon Compton subprocess gc — > 7c. The initial charm quark and the 
initial gluon are constituents of the initial hadrons. In addition, there is a leading order 
fragmentation process in which the photon is produced from quark or gluon fragmentation, 
e.g., gg cc followed by c — 7X, or qc qc followed by g ^ 7. At next-to- leading 
order in QCD, several subprocesses contribute to the ■y + c + X final state: gc gcj, 
gg — > CC7, qq cc'j, qc qcj, qc — > qcj, cc — > CC7, and cc — > CC7. A full next-to-leading 
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order calculation requires the computation of the hard-scattering matrix elements for these 
two-to-three particle production processes as well as the one-loop 0{al) corrections to the 
lowest order subprocess gc -yc. 

We are interested ultimately in the fully differential two-particle inclusive cross section, 
E^Ecda/d^p^d^Pc, where {E,p) represents the four- vector momentum of the 7 or c quark. 
For each contributing subprocess, this calculation requires integration over the momentum of 
the unobserved final parton in the two-to-three particle subprocesses {g, c, q, or q). Collinear 
singularities must be handled analytically by dimensional regularization and absorbed into 
parton momentum densities or fragmentation functions. In the theoretical analysis reported 
here, a combination of analytic and Monte Carlo integration methods is used to perform 
phase-space integrations over unobserved final-state partons and the momenta of the initial 
partons. This approach facilitates imposition of photon isolation restrictions and other se- 
lections of relevance in experiments. We work in the massless approximation, rric = 0. To 
warrant use of perturbation theory and the massless approximation, we limit our consider- 
ations to values of transverse momenta of the photon and charm quark pj,'^ > 5 GeV. 

In the lowest order direct subprocess, gc ■yc, the prompt photon emerges in isolation 
from the only other particle in the hard scattering, the charm quark. Long-distance quark- 
to-photon and gluon-to-photon fragmentation processes have been emphasized theoretically 
and parametrized phenomenologically in leading order |0] , and evolved in next-to- leading 
order PJ^. These terms may account for more than half of the calculated inclusive single 
photon cross section at modest values of transverse momentum at collider energies. Photons 
originating through fragmentation are likely to emerge in the neighborhood of associated 
hadrons. An experimental isolation restriction is needed before a clean identification can be 
made of the photon and a measurement made of its momentum. Isolation reduces the size 
of the observed fragmentation contribution. Photon isolation complicates the theoretical 
interpretation of results, however, since it threatens to upset the cancellation of infra-red 
divergences in perturbation theory In this paper, we include the fragmentation contri- 



butions, and we impose isolation requirements through our Monte Carlo method. 



A combination of analytic and Monte Carlo methods similar to that we employ in this 
paper has been used to carry out next-to-leading order QCD calculations of other processes 
including unpolarized and polarized inclusive prompt photon |TT| . [T^ , unpolarized and po- 
larized prompt photon plus jet production ||l3| , p!^ and unpolarized and polarized photon 



pair pr5||T^ production in hadron collisions, single and pair production of heavy gauge 



bosons ||T8l, and in our earlier work on spin-averaged 7c production The combination 
of analytic and Monte Carlo techniques used here to perform the phase space integrals is 
documented and described in detail elsewhere [|Tli|r5Hl7i|l8t ^|T6i with some details specific 



to the polarized case discussed in Ref. |T6|. We refer readers to those papers for further 
details. 

After theoretical expressions are derived in perturbative QCD that relate the spin- 
dependent cross section at the hadron level to spin-dependent partonic hard-scattering ma- 
trix elements and polarized parton densities, we must adopt models for spin-dependent par- 
ton densities in order to obtain illustrative numerical expectations. For the spin-dependent 
gluon density that we need, we use the three different parametrizations of AG{x,fi) sug- 
gested in Ref. . We generate the spin dependent polarized charm quark densities Ac(x, /i) 



perturbatively from these polarized gluon densities, beginning with the assumption that 
Ac(x, /io) = at the starting value for evolution /io = 1.5 GeV. Our polarized gluon and 
charm quark densities satisfy expected positivity constraints. 

We present spin-averaged cross sections differential in the transverse momenta and rapidi- 
ties of the photon and charm quark in pp collisions at the center-of-mass energy v^ = 200 
GeV typical of the Brookhaven RHIC collider. Our results on the longitudinal spin depen- 
dence are expressed in terms of the two-spin longitudinal asymmetry An, defined by 

^7,cc I I ^ _ (ji,c(j ) 

All = ' ' , (1.1) 

where -|-, — denote the helicities of the incoming protons. 

In Section II, we describe the next-to-leading order calculation of the spin- dependent 
cross section in perturbative QCD. Parametrizations of the spin-dependent gluon and charm 



quark densities are discussed in Sec. III. Differential cross sections and other numerical results 
are presented in Section IV. Summary remarks are collected in Section V. An appendix is 
included in which we derive analytic expressions for some of the parton level spin-dependent 
cross sections. 

II. CONTRIBUTIONS THROUGH NEXT-TO-LEADING ORDER 

The two-particle inclusive hadron reaction hi + h2 —>■ j + c+X proceeds through partonic 
hard-scattering processes involving initial-state light quarks q and gluons g. In lowest-order 
QCD, at 0{as), the only direct partonic subprocess is c -|- ^ 7 -|- c. In addition, there 
is a leading order fragmentation process in which the photon is produced from quark or 
gluon fragmentation, e.g., gg — > cc followed by c — > 7X, or qc — > qc followed by — > 7. 
Calculations of the cross section at order 0{a1) involve virtual gluon loop corrections to the 
0{as) direct subprocess as well as real gluon radiation contributions from a wide range of 
2 — > 3 parton subprocesses (of which some examples are shown in Fig. 1(c).) 

The full set of three-body final-state subprocesses is: 



g + c^ g + c + ^ (2.1a) 

g + g^c+c + -f (2.1b) 

q + q^c + c + ^ (2.1c) 

g + c^g + c + 7 (2. Id) 

q + c^q + c + ^ (2.1e) 

c + c^ c + c + -f (2. If) 

c + c^ c + c + ^ (2-lg) 



The physical cross section is obtained through the factorization theorem, 

~ i:/<*«.fe.A./)/^>..A./)^fep..y.* w). (2.2) 
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It depends on the hadronic center-of-mass energy 5", the transverse momenta and pj^ of 
the charm quark and photon, the rapidities yc and y^, and the relative azimuthal angle (p; 
/i/ is the factorization scale of the scattering process. The usual Mandelstam invariants in 
the partonic system are defined by s = (pi + P2Y, t = {pi — p-y)^, and u = {p2 — ^7)^, 
where pi and p2 are the momenta of the initial state partons and p^ is the momentum of the 
final photon. The indices ij label the initial parton channels whose contributions are added 
incoherently to yield the total physical cross section. The spin-averaged parton densities are 
denoted fh{x, fif). The partonic hard-scattering cross section a]f{s,pT, y, (p', fJ'f) is obtained 
commonly from fixed-order QCD calculations through 

^2^1,c r72fTT''='(^) r/2fTX'^'(^) 

The tree, virtual loop, and real emission contributions are labeled (a), (b), and (c). The 
parameter /i is the renormalization scale. 

The longitudinal spin-dependent cross section for hi + h2 — > 7 + c + X is very similar to 
Eq. (13): 

^ ''hih2 



dp^dp'^dycdy^dcj) 



/set lXct ■ ■ 
dxidx2Afl^ (xi,/i/)A/^^(x2,/i/) (g,PT, y, 0; /i/)- (2.4) 

The polarized parton densities are defined by 

Afiix, fif) = fl+{x, fif) - fl_{x, fif), (2.5) 

where flj.{x,fif) is the distribution of parton of type i with positive (+) or negative (— ) 
helicity in hadron h. Likewise, the polarized partonic cross section Aa'^'^ is defined by 

Aa^'^ = a^'^{+,+)-a^'%+,-), (2.6) 

with +, — denoting the helicities of the incoming partons. 



A. Leading order contributions 



In leading order in perturbative QCD, only one direct subprocess contributes to the 
hard-scattering cross section, the QCD Compton process eg —>■ 7c, unlike the case for 
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single inclusive prompt photon production, where the annihilation process qq 'yg also 
contributes. Since the leading order direct partonic subprocess has a two-body final state, 
the photon and c quark are produced with balancing transverse momenta. In addition, there 
are effectively leading-order contributions in which the photon is produced by fragmentation 
from a final-state parton. These are 

c + g ^ 9 + c 
g + g ^ c + c 
c + q c + q 
c + q ^ c + q 
c + c ^ c + c 
c + c —>■ c + c 

q + q^c + c. (2.7) 

If the photon is to be isolated from the observed charm quark, it arises from fragmentation 
of the gluon g and the non-charm quark q, respectively, in the cases of the first, third and 
fourth processes. In the other cases it is produced by fragmentation of one of the (anti) charm 
quarks. 

In a fully consistent next-to-leading order calculation, one should calculate the subpro- 
cesses in Eq. ( p.7|) to 0{al), since the photon fragmentation functions that are convoluted 
with the hard subprocess cross sections are of 0{aem/c(s)- For simplicity, we include them 
in 0(q;^) only. In fact, next-to-leading order fragmentation contributions to single prompt 



photon production have been included only once before p0| , p!2| . We expect the next-to- 
leading order corrections to the fragmentation contributions to be insignificant numerically 
especially after isolation cuts are imposed. This expectation is confirmed by a study per- 



formed in Ref. E^ on inclusive prompt photon production in the polarized case. It was also 



shown in Ref. |]2T| that predictions for the asymmetries are hardly affected by inclusion of 
the fragmentation contributions either before or after isolation of the photon. 
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B. Next-to-leading order contributions 



There are two classes of contributions in next-to-leading order. There are the virtual 
gluon exchange corrections to the lowest order process, eg 7c. Examples are shown in 
Fig. 1(b). These amplitudes interfere with the Born amplitudes and contribute at 0(aem«s)- 
They were calculated twice before in the spin-averaged case ||2^j23[] . At next-to- leading 
order there are also three-body final-state contributions, listed in Eq. (|2.1j ). Both the spin 
dependent and spin-averaged virtual loop and three-body matrix elements are taken from 
Ref. [^, where they are calculated for single inclusive prompt photon production. 



In Ref. p3| the original 'tHooft-Veltman-Breitenlohner-Maison (TVBM) p^ , p5| scheme 
was used to treat the 4-dimensional 75 matrix and antisymmetric e^jypo- tensor in n dimen- 
sions. These objects arise when the traces are taken to calculate the helicity dependent 
matrix elements. To project onto definite helicity states, h, for (anti-) quarks and A for 
gluons, the relations 

u{p,h)u{p,h) = ^Y{pKl-h-f5) (2.8) 

and 

1 r n^Ti''^" 

(2.9) 



p.p' 



are used. In the TVBM scheme, n-dimensional Minkowski space is divided into a 4- 
dimensional and a (n — 4)-dimensional part. Any vector k has a 4-dimensional, k, and 
a (n — 4)-dimensional part, k. This means that the n dimensional matrix elements calcu- 



lated in Ref. ||2^ contain products of these (n — 4)-dimensional 'hat' momenta. In Ref. ||23 
it was shown that these contributions are non-zero in the collinear limit only, and thus in 
our case they are included in the hard collinear integrals given in the appendix. A detailed 
account of all aspects of the phase space integration will be given elsewhere pH]. It should 



be mentioned here that the matrix elements were also calculated in Ref. in a scheme 
different from the TVBM scheme outlined here. Their results seem similar to those of Ref. 
23| , but a detailed comparison of the matrix elements has not yet been made. 
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The main task of our calculation is to integrate the three-body matrix elements over the 
phase space of the unobserved particle in the final state. The situation here is different from 
the standard case of single inclusive particle production because we wish to retain as much 
control as possible over the kinematic variables of a second particle in the final state, while at 
the same time integrating over enough of the phase space to ensure cancellation of all infrared 
and collinear divergences, inherent when massless particles are assumed. All the processes 
of Eq. (|2.1|) , except the first, involve collinear singularities but no soft singularities. The 
soft and collinear singularities are first exposed by integrating the three-body phase space 
in n-dimensions over the soft and collinear regions of phase space. These soft and collinear 
regions are separated from the rest of the three body phase space by introducing soft and 
collinear cut-off parameters 6s and Sc- The collinear and soft singularities are cancelled and 
factored in the MS scheme as explained in Ref. 0. 

At 0{al) there are, in addition, fragmentation processes in which the hard-scattering 
two-particle final-state subprocesses 

c + 5' ^ 7 + c 
c + c ^ 'y + g 



are followed by fragmentation processes c cX, in the case of the first subprocess, and 
g —>■ cX in the cases of the last two. These should be included because we have factored 
the collinear singularities in the corresponding three-body final-state processes into non- 
perturbative fragmentation functions for production of a charm quark from a particular 
parton. As a first approximation, we estimate these fragmentation functions by 



q + q-^l + g 



(2.10) 



2n 




(2.11) 



and 




(2.12) 
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where Pij{z) are the lowest order sphtting functions for parton j into parton i |2^; and 
asifi"^) is the strong couphng strength. 

For completeness we have gathered some of the main ingredients used in calculating the 
polarized cross section in the appendix. Similar results for the unpolarized case can be found 
in Ref. ||^. Further details and a complete list of the 4-dimensional three-body final-state 
matrix elements will be given in Ref. . 



III. POLARIZED PARTON DENSITIES 

Our spin-dependent gluon density and our spin- dependent up, down, and strange quark 
densities (and their antiquark components) are adapted from the set of polarized parton 
densities published by Gehrmann and Stirling (GS) [1^. However, this set of spin-dependent 



densities, as well as other recent parametrizations available for general use, does not include 
a charm contribution At the values of squared four- momentum transfer, Q^, at which 



spin-dependent deep-inelastic scattering data have been obtained, the charm contribution is 
expected to be very small. Fits to the data are thus insensitive to the charm contribution. 
The situation is somewhat different in the spin-averaged case, and mechanisms for inclusion 
of charm are in general use. 

In particular, a charm component is included in the CTEQ4M spin-averaged parametriza- 
tion [0 at a scale fio = itLc, where ~ 1.5 GeV is the charm quark mass. At this threshold 
scale the charm quark density is assumed to have zero value, c{x,fio) = 0. Non-zero charm 
quark and antiquark densities are generated at larger values of fi via the DGLAP pertur- 
bative QCD evolution equations through gluon splitting into cc pairs. Since we use the 



CTEQ4M distributions ||3T| as our spin-averaged parton densities, it is reasonable to use 
a similar mechanism to generate the spin- dependent charm densities, with spin-dependent 
Altarelli-Parisi splitting functions used in the evolution. 

The evolution mechanism outlined above implies that the spin-dependent charm quark 
density will depend strongly on the assumed value of the spin- dependent gluon density. 
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Uncertainty in knowledge of the gluon density will propagate to the charm density. The 
current deep inelastic scattering data do not constrain the polarized gluon density very 
tightly, and most groups present more than one plausible parametrization. Gehrmann and 
Stirling present three such parametrizations, labelled GSA, GSB, and GSC. In our 
parametrizations of the spin-dependent densities we begin with the GS parametrizations as 
our initial set at a scale fio = ttLc and we evolve them to obtain spin-dependent gluon and 
quark densities at greater values of /z, including a spin-dependent charm quark density. We 
use next-to-leading order DGLAP evolution equations with the number of flavors, Uf = 4. 



Details of the x-space evolution will be reported elsewhere Ref. ||32|. Our approach is to 
adopt the GS polarized parton densities as our starting values and to extend them to include 
charm. 

In Fig. 2(a), we display the ratio of our spin-dependent and spin-independent gluon 
densities, AG{x, fi)/G{x, fi), at the hard scale /i^ = 100 GeV^ for the GSA, GSB, and GSC 
choices. In the GSA and GSB sets, AG{x, fio) is positive for all x, whereas in the GSC set 
AG{x,fio) changes sign. After evolution to /i^ = 100 GeV^, AG{x,^) remains positive for 
essentially all x in all three sets, but its magnitude is small in the GSB and GSC sets. In 
Fig. 2(b), we show the spin-dependent charm quark distributions. Positivity is well satisfied, 
as it is for AG{x, fi)/G{x, fi). The value of Ac(x, yu)/c(x, /i) is reasonably substantial in the 
GSA case, reflecting the large size of AG{x, fi) from which Ac{x, fi) is generated. The size of 
Ac{x, fi)/c{x, fi) is correspondingly smaller in the GSB and GSC cases, hovering near zero 
in the GSC case. 

As mentioned, our charm quark densities depend on the choice made for the gluon den- 
sities. In addition, we assume that there is no intrinsic charm quark density. If an intrinsic 
density were adopted, one might begin with non-zero charm and anti-charm densities at fio 
whose x-dependences need not be assumed identical, in both the spin-averaged and spin- 
dependent cases. In a sense, our assumed forms for the spin-dependent charm quark densities 
represent the most conservative possibilities, and one should be alert experimentally to more 
interesting outcomes. 
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IV. NUMERICAL RESULTS 



In this section we present and discuss spin-averaged differential cross sections and two- 
spin longitudinal asymmetries for the joint production of a charm quark and a photon at 
large values of transverse momentum. All results are displayed for pp collisions at the 
center-of-mass energy = 200 GeV typical of the Brookhaven RHIC collider. To obtain 
the spin-averaged differential cross sections presented in this paper, we convolute our hard- 
scattering matrix elements with the CTEQ4M parton densities We use the standard 
two-loop formula for the strong coupling strength with four massless flavors of quarks. We 
set Ag^o = 0.202 GeV (the CTEQ4M value). The spin-dependent parton densities are 
described in Sec. III. 

Very similar differential distributions may be obtained if other parton sets are used 
instead, with quantitative differences reflecting differences among gluon and charm quark 
densities in the different sets 0]. We set the renormalization, factorization, and fragmen- 
tation scales to a common value fi = pj^ in most of our calculations. Since there are two 
observed particles in the final state, the charm quark and the photon, both of whose trans- 
verse momenta are large, an alternative choice might be /i = or some function of p^ 
and p^. The results of our calculations show that the magnitudes of pj^ and p^ tend to 
be comparable and that dependence of the cross sections on /i is slight. Therefore, choices 
of /i different from fi = p^ should not produce significantly different answers, and we have 
verified this supposition in representative cases. 

In collider experiments a photon is observed and its momentum is well measured only 
when the photon is isolated from neighboring hadrons. In our calculation, we impose isola- 
tion in terms of the cone variable R: 

.^(Ay)2 + (A0)2 < R. (4.1) 

In Eq. ( [4.1|) , Ay (A0) is the difference between the rapidity (azimuthal angle in the transverse 
plane) of the photon and that of any parton in the final state. The photon is said to be 
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isolated in a cone of size R if the ratio of the hadronic energy in the cone and the transverse 
momentum of the photon does not exceed e = 2GeV/pj.. We show distributions for the 
choice -R = 0.7. 

In Fig. 3(a) we show the spin- averaged differential cross section as a function of the 
transverse momentum of the photon pj,, having restricted the transverse momentum of the 
charm quark to the range > 5 GeV. The rapidities of the charm quark and photon are 
restricted to the central region, — 1 < y'^''^ < 1 in order to mimic the central region coverage 
of the major detectors at RHIC. The solid curve shows our prediction when contributions 
are included from all subprocesses through 0{a'^). The dominance of the eg subprocess is 
illustrated by our dashed curve. This dominance is the basis for the statement that the 
spin-averaged cross section at collider energies can be used to determine the magnitude and 
Bjorken x dependence of the spin- averaged charm quark density 0,0. 

In Figs. 3(b) and (c), we present the two-spin longitudinal asymmetries. All, as a func- 
tion of pj', for the same kinematic selections as made for Fig. 3(a). Results are shown for 
three choices of the polarized gluon density (and, correspondingly, for the polarized charm 
quark density). The asymmetry becomes sizeable for large enough only in the case of 
the GSA parton set. 

As noted above the eg subprocess dominates the spin-averged cross section. It is in- 
teresting and important to inquire whether this dominance persists in the spin-dependent 
situation. The dot-dashed curve in Fig. 3(b) shows the contribution to the asymmetry from 
the polarized eg subprocesses in the case of the GSA set. It is positive for all pj^. At small 
Pj', the net asymmetry is driven negative by a large contribution from the gg subprocess. For 
this GSA set, we see that once it becomes sizable (e.g., 5% or more), the total asymmetry 
from all subprocesses is dominated by the large contribution from the eg subprocess. In 
Fig. 3(c), using a more expanded scale, we replot the overall asymmetry for the GSC choice 
of the polarized gluon density, and we also show the contribution from the eg subprocess 
alone. In this case, the overall asymmetry All itself is small, and the contribution from the 
eg subprocess cannot be said to dominate the answer. 
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In interpreting the results presented in Figs. 3(b) and (c), we begin with the supposition 
that the polarized gluon density will have been determined from data on inclusive prompt 
photon production. The question to pose is whether asymmetries of the type shown in 
Figs. 3(b) and (c) could shed light on the polarization of the charm quark density. If a 
large asymmetry is measured, similar to that expected in the GSA case at the larger values 
of p^, Fig. 3(b) shows that the answer is dominated by the eg contribution, and the data 
will serve to constrain Ac(x, /i). On the other hand, if AG(x, /i) is small and a small 
asymmetry is measured, such as shown in Fig. 3(c), one will not be able to conclude which 
of the subprocesses is principally responsible, and no information could be adduced about 
Ac(x, fi). 

In Fig. 4 we show distributions in the rapidity of the charm quark, y'^, for —1 < y'^ < 1, 
> 5 GeV, and 4 GeV < < 50 GeV. The spin-averaged distribution in y'^, shown in 
Fig. 4(a) is fairly broad, with full-width at half-maximum of nearly 3 units in rapidity. The 
dashed curve in Fig. 4(a) shows that the contribution from the eg subprocess is dominant. 
The asymmetries are shown in Fig. 4(b) as functions of y'^. They are small for all three 
choices of the the polarized gluon density. These results are consistent with those shown in 
Fig. 3(b), although they might not appear to be at first glance. One must bear in mind that 
it is the region of relatively small pj^ that dominates the integral over pj,. The asymmetries 
shown in Fig. 3(b) are small at modest values of pj. In the GSA case, the total asymmetry 
at relatively small values of pj, is opposite in sign to the positive contribution from the eg 
subprocess. 

The structure of the QCD hard-scattering matrix element produces positive correlations 
in rapidity at collider energies. To examine correlations more precisely, we study the 
spin-averaged cross section and the asymmetry as functions of the difference of the rapidities 
of the photon and charm quark. Results are shown in Fig. 5. 
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V. SUMMARY AND DISCUSSION 



In this paper we present the results of a calculation of the longitudinal spin-dependence of 
the inclusive production of a prompt photon in association with a charm quark at large values 
of transverse momentum. This analysis is done at next-to-leading order in perturbative 
QCD. We employ a combination of analytic and Monte Carlo integration methods in which 
infrared and coUinear singularities of the next-to-leading order matrix elements are handled 
properly. We provide differential cross sections and polarization asymmetries as functions of 
transverse momenta and rapidity, including photon isolation restrictions, that may be useful 
for estimating the feasibility of measurements of the spin-averaged and spin-dependent cross 
sections in future experiments at RHIC collider energies. We show that the study of two- 
particle inclusive spin-averaged distributions, with specification of the momentum variables 
of both the final prompt photon and the final heavy quark, tests correlations inherent in 
the QCD matrix elements and should provide a means for measuring the charm quark 
density in the nucleon [Q. In the spin-dependent case, significant values of An (i.e., greater 
than 5 %) may be expected for pj, > 15 GeV if the polarized gluon density AG{x,fi) is 
as large as that in the GSA set of polarized parton densities. If so, the data on associated 
production could be used to determine the polarization of the charm quark density in the 
nucleon. On the other hand, for small AG{x, /i), dominance of the eg subprocess is lost, and 
Ac(x, /i) is inaccessible. 

Our spin-averaged and polarized charm quark densities are generated perturbatively 
from gluon splitting into charm quark and charm antiquark pairs. They depend entirely on 
the choice made for the gluon densities: a small polarized gluon density leads to a small 
polarized charm density. Reality may well be different. For example, there could be a 
non-perturbative intrinsic charm quark component. In a sense, our assumed forms for the 
spin-averaged and spin-dependent charm quark densities represent the most conservative 
possibihties. 

In a typical experiment, the momentum of the quark may be inferred from the momen- 
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turn of prompt lepton decay products or the momentum of charm mesons, such as D*'s. 
Alternatively, our distributions in may be convoluted with charm quark fragmentation 
functions, deduced from, e.g., e'^e~ annihilation data, to provide distributions for the prompt 
leptons or D*'s. 
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APPENDIX A: ANALYTIC EXPRESSIONS 



In order to make this paper reasonably self-contained, we collect in this appendix all the 
formulae we use in the calculation of the polarized cross section. We label the momenta for 
the generic three-body final-state process by 

P1+P2 ^P3+P4+P5, (Al) 

where ps denotes the photon and denotes the observed charm quark, and we define the 
usual Mandelstam invariants 

s = {pi+P2f. (A2) 
We express the two-body final-state cross sections in terms of the scaled variable v, where 

v^l + i. (A3) 



1. Two-body contributions 

The effective two-body contribution includes the 0{a'^) virtual gluon-exchange loop con- 
tributions and the soft and/or coUinear parts of the three-body contributions (this remark 
applies to initial-state coUinear contributions only, as explained later) . After all soft pole sin- 
gularities are cancelled and all coUinear pole singularities are factored, the effective two-body 
contribution is expressed as 



dvdxidx2 



dv 



dAa 



HO 



dv 



-{eg 7c) 



(A4) 



plus terms in which the beam and target are interchanged. We use subscript c to label the 
charm (or anti-charm) quark. 
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We define 



AT, 



eg 



(A5) 



and Vi = 1 — V. In Eq. 

■{eg 7c) 



dv 



sNc V 
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[-UCpAT^g + 2{2Cf + iVj In^ 5,AT,, 



-A^^ In ^ATc„ + Cf(3 + 41n5J ( In ^ + In ^ | AT,„ 



iV ox 

+ + 12 In 5,) In ^AT,, + -(llA^, - 2Ar^) In —AT,, 



+ {2Cf - A'c) In^ ^^ATcg + 4 In 5s{N^ In i;i + 2Cf - In t;) AT, 



eg 



+ 2(2Cf - Nc) In vi In t;ATcg - 2(:7^7r 



2{3-4:V-v') 



3fi 



+ A^cvr- 



, (3 - 2w - 2v^ 
3vi 



+ In^ NJl + v)- 2Cf 



+ \Yivi\ 2Nc + 2Cf^^-^^] + 2{2Cf - N,)AT^gU2{vi) 



Vl 



+ 2Ar,AT,3Li2(t;) 



(A6) 



The scales M and M" are the factorization and fragmentation scales, respectively, on the 
initial parton and final-state charm quark legs, and fi is the renormalization scale. Ci? = 4/3 
is the quark-gluon vertex color factor, A'c = 3 is the number of colors, Cc is the fractional 
charge of the charm quark, 6s and 6c are the soft and collinear cutoff parameters introduced in 
Sec. II, Li2(x) is the dilogarithm function, and Oem is the electromagnetic coupling constant. 

The remnants of the factorization of the hard collinear singularities are 



dv 



2s N, 



c 



eg 



i-<5= dz 



fX2 



X2 



-A/f (^^,M^)AP,,(z) 



+ 



dz 



+ 



-A/, 



-A/-rA,M^ AP,,(. 



A I 



1 APggiZ 



xi Z 



(AT) 



The last distribution,A/g^(xi, M^), in Eq. impl ies a sum over the flavors of quarks from 



the cc, eg, and cq initial states. The remaining two processes, cc and qq, do not have initial 
state coUinear singularities and thus do not contribute to this part of the cross section. 
The polarized splitting functions APjj, are 

'1- z 



-6,— ]-APr{z). 



(A8) 



z "M2^ 

The functions APij{z,e), the spin dependent Altarelli-Parisi splitting functions in 4 — 2e 
dimensions, are 

'1 + z^ 



AP,,{z,e) = CF 



1-z 



+ 3e(l - z) 



AP,,iz,e)^-[i2z-l)-2eil-z)] 



AP,,{z,e) = 2Nc 



2z + l + 2e(l-z] 



1-z 

APgg{z, e) = Cf [(2 -z)+ 2e(l - z)] . 
The functions APlj{z) are defined by the relation 



(A9) 



AP,j{z,e)^APij{z)+eAPUz). 



(AlO) 



2. Pseudo- Two-Body Contributions 

Since we are interested in distributions in the kinematic variables of two final-state 
partons, the photon and the c quark, we can define variables that depend on the momenta 
of both. An example is the variable z, defined by 

- = -f#- (All) 
\Pt\ 

Whenever there is a third parton in the final state, the distribution in z (or in other analogous 
variables) will differ from a 6 function when the third parton carries a finite momentum, 
even if it is coUinear to one of the other final partons. For this reason we designate as 
"pseudo-two-body contributions" those for which the third parton is coUinear to either the 
final photon or the charm quark. These contributions are expressed, respectively, by the 
equations 
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-{ah cq)dxidx2dzdv, 



abcq ' 

dAa 



2-n 



M'2 



dv 



(A12) 



and 



Act c/co// 



ahd 



112 



dAa 

dv 



{ah — > 'yd)dxidx2dzdv. (A13) 



The functions P.y/q{z) and P^^g{z) are the quark-to-photon sphtting function and 0(e) piece, 
respectively. They have the same form as Pgg, with the color factor replaced by the square 
of the quark charge. The scale M' is the fragmentation scale for quark fragmentation into 
a photon. In Eq. (|A13|) , 



PUz,M"^)=PUz) In 



z{l — z)6c 



M"2 



(A14) 



where Pcd{z) represents the splitting functions Pqg{z) and Pqq{z). The unpolarized splitting 
functions are 



Pqq{z, e) = Cf 



l + z' 



-e{l-z) 



Pgg{z,e)=2Nc 



1 - z 

z' + {l- zY - e 
1 - z 



Pgq{z, e) — Cp 



1 — z z 
l + {l-zf 



+ z{l-z) 



ez 



(A15) 



The functions PL{z) are defined by the relation 



P,,{z,e) = P,,{z)+ePr{z). 



(A16) 



In Eq. ( |A12|) , q can be a charm quark or a charm anti-quark, or a(an) (anti)quark of any 
flavor in the case of cq cq. In Eq. (|A13|) , d is either a gluon, a charm quark, or a charm 
antiquark. Note that the splitting functions P-y/q and Pcd are the usual spin averaged and 
not the spin dependent ones, since the final state is not polarized. 
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These contributions are referred to usually as the remnants of the factorization of the 
hard collinear singularities and are regarded as two-body processes, or as parts of the frag- 
mentation contributions because of their dependence on the factorization scales. We prefer 
to regard them as pseudo-two-body contributions. When we examine either the charm quark 
or photon momentum distributions, these contributions populate the same regions of phase 
space as the other three-body contributions in Eq. (|A20|) , unlike the effective two-body con- 
tributions. The pseudo-two-body contributions are usually negative in overall sign due to 
the large logarithms of the cut-off parameters, as are the two-body contributions discussed 
above. 



3. Photon Fragmentation Contributions 

As mentioned in Sec. II, we include the quark-to-photon and gluon-to-photon fragmen- 
tation contributions at leading order only. We convolute the 2 — 2 hard scattering subpro- 
cess cross sections for the processes listed in Eq. ( p.7|) with photon fragmentation functions 
D^li{z, M'^); M' is the fragmentation scale, the same scale at which we subtract the collinear 
singularities on the photon leg of the three-body processes. The expression for the cross sec- 
tion is 

abj 

The matrix elements for the hard subprocess cross section can be found, for example, in 
Ref. ra. 



^^frag/^ = J ^fai^u M^)Af,''{x2, M^)D^/j{z, M'^)^{ab J c)dx,dx2dzdv . (Al7) 



4. Charm Fragmentation Contributions 

In integrating some of the three-body matrix elements over phase space we encounter 
configurations in which the produced charm quark is collinear with an anti-charm quark or 
a gluon in the final state. These situations lead to a collinear singularity in the massless 
approximation. They occur for the processes of Eq. ( p.laj) , (p.lc|) , and (p.lf|) . We factor 
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these singularities into a fragmentation function Dc/i{z, M"^) for parton i to produce a charm 
quark with momentum fraction z. The contributing subprocess cross sections are 



dAa . _ . 2Cf TiaasC^ f v I -v 

—i—m 19) = — — + 

dv Nc s V 1 — f 

dAa. . naascl f 1 — (1 — v)"^ 



A^Wc = E / Af^{xuM^)Af,''{x2,M^)D,/a{z,M"^)^{ab ^ ^d)dx,dx2dzdv. (A19) 

abd 



la) = 


2Cf 


Nc 


-fq) = 




Ncs 


;iven by 





(A18) 



5. Three-body contributions 

The non-colhnear three-body final-state contributions are calculated from the expression 
Aa^body = Y.j ^fa i^i, M2)A/,^(x2, M^)dAaiab ^ ^cd)dxidx2dn, (A20) 

abd 

with Q representing the angles and other variables that are integrated over. Whenever an 
invariant Sij or tij falls into a collinear or soft region of phase space, that contribution from 
the subprocess is excluded. The three-body contribution shows no dependence on the fac- 
torization scale of the final-state charm or photon legs, although we have factored collinear 
singularities at scales M" and M', respectively, on these legs of the three-body subprocesses. 
However, Eq. ( |A20| ) does contain implicit logarithmic dependence on the soft and, in partic- 
ular, the collinear cutoffs discussed in Sec. II. Both collinear cutoff and factorization scale 
dependences are contained in the pseudo-two-body contributions discussed above. 
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FIGURE CAPTIONS 



[1] (a) Lowest order Feynman diagrams for 7 plus c quark production; ki and k2 
are the four- vector momenta of the photon and charm quark, (b) Examples of 
virtual corrections to the lowest order diagrams, (c) Examples of next-to-leading 
order three-body final-state diagrams for the gc initial state. 

[2] (a) Ratio of the polarized to the unpolarized gluon distribution AG{x, //) / G{x, //) 
at is^ = 100 GeV^ for the GSA, GSB and GSC spin-dependent parton dis- 
tributions and the CTEQ4M spin-averaged distributions. Evolution begins at 
fjLo — 1-5 GeV, and the number of flavors n/ is set to 4 in the evolution, (b) 
Same as (a) but for the charm distributions. The charm quark density is gener- 
ated from bremsstrahlung of gluons into charm-anticharm pairs. 

[3] (a) Spin-averaged cross section da/dpj^ as a function of the transverse momen- 
tum of the photon for p -|- p ^ ^ + c + X at = 200 GeV. The transverse 
momentum of the charm quark is restricted to pj^ > 5 GeV, and the rapidities 
of the photon and charm quark are restricted to the interval —1.0 < y < 1.0. 
Two curves are drawn. The solid curve shows the cross section that results from 
the inclusion of all subprocesses through 0(q;^), and the dashed curve shows the 
contribution from the dominant eg subprocess, only, through O(a^). (b) For 
the same selections as in (a), the two-spin longitudinal asymmetry is shown for 
the three choices of the spin-dependent gluon density. In the case of GSA, the 
dot-dashed curve shows the contribution to A^^^ from the eg subprocess only, (c) 
For the same selections as in (a), the two-spin longitudinal asymmetry All is 
shown for the GSC spin-dependent gluon density (dotted hue), along with the 
contribution from the eg subprocess only. 

[4] (a) Spin-averaged cross section da/dy"^ as a function of the rapidity of the charm 
quark f or p -|- p — > 7 -|- c -|- X at = 200 GeV. The photon rapidity is restricted 
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to —1.0 < y'^ < 1.0, and 4 < < 50 GeV. The transverse momentum of the 
charm quark is selected to be in the region > 5 GeV. The sohd curve shows the 
cross section that results from the inclusion of all subprocesses through 0{al), 
and the dashed curve shows the contribution from the dominant eg subprocess 
through 0(q;^). (b) For the same selections as in (a), the two-spin longitudinal 
asymmetry is shown for the three choices of the spin-dependent gluon density. 
The sohd, dashed, and dotted curves represent the results from calculations based 
on our modified GSA, GSB, and GSC densities, respectively. In the GSA case, 
the dot-dashed curve shows the contribution to All from the eg subprocess only. 

[5] (a) Spin- averaged cross section da/dAy as a function of the difference Ay — 
y'^ — y^ of the rapidities of the photon and charm quark, foicp + p^^ + e + X 
at = 200 GeV. The transverse momentum of the photon is selected to be in 
the interval, 4 < pj, < 50 GeV, and its rapidity is restricted to —1.0 < y"' < 1.0. 
The transverse momentum of the charm quark satisfies > 5 GeV. The solid 
curve shows the cross section that results from the inclusion of all subprocesses, 
and the dashed curve shows the contribution from the dominant eg subprocess 
only, (b) For the same selections as in (a), the two-spin longitudinal asymmetry 
is shown for the three choices of the spin-dependent gluon density. The sohd, 
dashed, and dotted curves represent the results from calculations based on our 
modified GSA, GSB, and GSC densities, respectively. 
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